Nonlinear theory of fractional microwave-induced magnetoresistance oscillations 
in a dc-driven two-dimensional electron system 
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Microwave-induced nonlinear magnetoresistance in a dc-driven two-dimensional electron system 
is examined using a multi- photon-assisted transport scheme direct controlled by the current. It is 
shown that near the 2nd subharmonic of the cyclotron resonance, the frequency of the resistivity 
oscillation with the magnetic-field-normalized current-density is double that at the cyclotron reso- 
nance and its harmonics, in excellent agreement with recent experimental findings by Hatke et al. 
[Phys. Rev. Lett. 101, 246811 (2008)]. The current-induced alternative emergence of resonant 
two-photon and single-photon processes is responsible for this frequency doubling. Near the third 
subharmonic of the cyclotron resonance, the current-induced consecutive appearance of resonant 
0-/3- photon, two-photon, and single-photon processes may lead to the frequency tripling of the 
resistivity oscillation. 

PACS numbers: 73.50.Jt, 73.40.-c, 73.43.Qt, 71.70.Di 
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I. INTRODUCTION 

A number of extraordinary magnetotransport phe- 
nomena of irradiated two-dimensional (2D) elec- 
trons in very high Landau levels (LLs), especially 
the microwave- induced resistance oscillation (MIRO) 
and the zero-resistance state, has been the focus 
of intensive experimentall^^'^'^i'^'a'ia^i^^^ and 
theoreticali^^ii^ia^i^i^^^^^^S^^ studies 
in the past few years. 

Under the irradiation of a moderate microwave of fre- 
quency w, the most striking resistance oscillations show 
up in the linear photoresistivity as a function of the in- 
verse magnetic field 1/-B, featuring the periodical ap- 
pearance of peak- valley pairs around uj/lOc — 1,2,3... 
{lOc = eB/m is the cyclotron frequency and m the elec- 
tron effective mass). The cyclotron resonance and its 
harmonics are node points of peak- valley pairs, where the 
photoresistivity vanishes. In addition, secondary peak- 
valley pairs were also observed in the linear resistance 
near certain fractional values, lo/lUc — 1/2, 1/3,2/3, and 
3/2. These fractional microwave-induced resistance os- 
cillations were referred to real multiphoton participating 
processesi^ They were also referred to higher harmonic 
effect)^ the resonant series of consecutive single-photon 
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microwave-induced sidebands De- 



spite different scenarios have been proposed for these 
fractional MIROs of linear photoresistivity, they have so 
far been theoretically demonstrated to clearly appear in 
a real system only within the multiphoton-assisted scat- 
tering model^i. 

Another notable effect recently discovered is the signif- 
icant role of a dc current. A relatively weak current alone, 
without microwave radiation, can also induce substantial 
magnetoresistance oscillations and zero-differential resis- 
tance states ,^^i22i^i2adaiiiii^d^ Simultaneous applica- 
tion of a finite dc current and a microwave radiation leads 
to very interesting and complicated oscillatory behavior 
of resistance and differential resistance4^i^i^i^i^ These 



facts clearly indicate that microwave-irradiated nonlin- 
ear magnetotransport is remarkably different from the 
linear one. So far experimental and theoretical investiga- 
tions on nonlinear MIROs have been carried out in the 
magnetic field range lower than the cyclotron resonance 
and never anticipated any anomaly elsewhere. Very re- 
cent experiments by Hatke et ali^ disclosed sharply dis- 
tinct nonlinear transport behavior in the subharmonic 
fields of the cyclotron resonance. Particularly, they found 
that at the second cyclotron resonance subharmonic, 
uj/lUc = 1/2, the differential resistivity oscillates as a 
function of the current density in the same way as that 
at the cyclotron resonance, ui/ujc — 1/2, despite a fac- 
tor of two difference in magnetic fields. In other words, 
at u/lOc = 1/2 the frequency of the resistivity oscilla- 
tion with the magnetic-field-normalized current-density 
is twice that at the cyclotron resonance and its harmon- 
ics. This is in vast contract with the observed behavior 
of the resistivity oscillations in the magnetic field regimes 
of cyclotron resonance and its harmonics (uj / lOc > 1) ■ A 
nonlinear theoretical analysis capable of identifying the 
specific mechanism responsible for fractional MIROs and 
frequency doubling is highly desirable. 

This paper presents such an analysis on microwave- 
induced nonlinear magnetotransport in high Landau lev- 
els in GaAs-based 2D heterosystems. With the micro- 
scopic balance-equation scheme direct controlled by the 
current, we are able to constitute a unified description 
for both the cyclotron resonance harmonic and subhar- 
monic regimes of the magnetic fields and predict the re- 
sistivity oscillations with applied current density in ex- 
cellent agreement with the experimental observation in 
both regimes. 



II. FORMULATION 

Our examination is based on a current-controlled 
scheme of photon-assisted transport, which deals with a 
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2D system of short thermalization time having Ng elec- 
trons in a unit area of the x-y plane. These electrons, 
subjected to a uniform magnetic field B = (0, 0, B) 
in the z direction and scattered by random impurities 
and by phonons in the lattice, perform an integrative 
drift motion when an incident microwave electric field 
EisSmujt irradiates on and a current flows in the plane. 
In terms of the center-of-mass (cm.) momentum and 
coordinate defined as P = J^j Pj\\ ^^'^ ^ = 12j '"i 
with PjIi = {pjx,Pjy) and rj = {xj,yj) being the mo- 
mentum and coordinate of the jth electron in the 2D 
plane, and the relative electron momentum and coordi- 
nate p^.ji = Pj Ii —P/Ns and r'^ = rj —R, the Hamiltonian 
of this system can be written as the sum of a cm. part 
Hem. and a relative electron part -fferr^ 

Hem. = TTT^ (P - N,eA{R)y ~ N.eiEo + Et) ■ R,{1) 
+ Y,VM-r',,z,,z,), (2) 

Kj 

together with electron-impurity and electron-phonon 
couplings Hci and i/cp- Here A{r) = {—By, 0) is the 
in-plane component of the vector potential of the per- 
pendicular magnetic field, Eq and Et are the dc and ac 
components of the uniform electric field inside, m and 
are, respectively, the electron effective mass parallel and 
perpendicular to the plane, V{z) and Vc(r- — r^, Zi,Zj) 
stand for the confined and Coulomb potentials. 

The separation of the electron Hamiltonian into a cm. 
part and a relative electron part amounts to look at elec- 
trons in the reference frame moving with their center of 
mass. The most important feature of this separation is 
that a spatially uniform electric field shows up only in 
Hem., and that Ha- is the Hamiltonian of a many parti- 
cle system subject to a magnetic field without the electric 
field. This enables to deal with the relative electrons in 
the magnetic field without tilting the Landau levels. 

We proceed with the Heisenberg operator equations for 
the rate of change in the cm. velocity V = —i[R,H]: 
V = —i[V,H], and that of the relative electron energy 
Her'. Her — —i[Hcr,H]. The cm. coordinate R and 
velocity V in these equations can be treated classically, 
i.e., as the time-dependent expectation values of the cm. 
coordinate and velocity, R{t) and V{t), such that R{t) — 
R{t') = Jj* V{s)ds. In the case of the steady transport 
under a modest radiation of single frequency uj it suffices 
to assume that the cm. velocity, i.e., the electron drift 
velocity, consists of a dc part v and a stationary time- 
dependent part of the form 

V{t) = V + Vi cos((jjt) + V2 sm{ujt). (3) 

For high mobility and high carrier density 2D systems 
at low temperatures where effects of electron-impurity 



and electron-phonon scatterings are weak in comparison 
with internal thermalization of this many electron sys- 
tem, it is a good approximation to carry out the statis- 
tical average of the above operator equations to leading 
orders in Hd and H^p. For this purpose we only need 
to know the distribution of relative electrons without be- 
ing perturbed by Hd or H^p. The distribution function 
of the system described by Hamiltonian ^ without elec- 
tric field, is an isotropic Fermi-type with a single temper- 
ature Tg. Such a statistical average of the above operator 
equations yields the following force and energy balance 
equations in the steady state: 

NscEo + Nse{v X B) + F = 0, (4) 
N.eEo -v + Sp-W ^0. (5) 

In this, 

oc 

^-E|f^(9||)r E Q||J'(0^f2(9||,^o-nu;) (6) 

q|l n— — oo 

is the time-averaged damping force against CM motion, 
Sp is the time-averaged rate of the electron energy- 
gain from the ac field, having an expression obtained 
from Eq. (6) by replacing gn factor by nui. In Eq. ©, 
C/(q||) is the effective impurity potential, JniO the 
Bessel function of order n, n2{q\\,fl) is the imaginary 
part of the electron density correlation function at elec- 
tron temperature in the presence of the magnetic 
field, and ujo = Q\\ ■ i>. The argument in J„(^) is 
£, = y^iq\\ ■ vi)'^ + (q|| • v-i)^ jijj. Note that, although con- 
tributions of phonon scattering to F and Sp are neglected 
here, it provides the main channel for electron energy dis- 
sipation to the lattice with a time-averaged energy-loss 
rate W , having the expression as given in Ref.fTol. 

The ac components v\ and V2 of the electron drift ve- 
locity should be determined selfconsistently from the in- 
cident ac field E\^ by the Maxwell equations connect- 
ing both sides of the 2D system, taking account the 
scattering-related damping forces F^ and Fc^^ However, 
for the high-mobility systems at low temperatures, the 
effects of these scattering-related forces are much weaker 
than that of radiative decay^ and always negligible, 
hence V\ and Vi are directly determined by the setup 
of the 2D system in the sample substrateJ^ 

The effect of interparticle Coulomb interaction is in- 
cluded in the density correlation function to the degree 
of electron level broadening and screening (considered in 
the effective impurity and phonon potentials). The re- 
maining 1/2(911, ^) function in Eqs. (O is that of a non- 
interacting 2D electron gas in the magnetic field, which 
can be written in the Landau representation as^^ 

772(911, /2) = ^ E C„,„'(?|(Zf|/2)772(n, n', Q), (7) 

ir2(n,n',/2) = -^ y d£ [/(£)-/(£ + /?)] 

X ImG„(£-l- /2)ImG„'(e), (8) 
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where = y/T/\eB\ is the magnetic length, C„,„+i(F) = 
nl[{n + l)\]-^Y^e-^[Ll{Y)]^ with Ll{Y) the ' associate 
Laguerre polynomial, f{e) = {exp[(e — A^)/7e] + 1}^^ 
is the Fermi function at electron temperatm'c Tg, and 
ImG„(£) is the density-of-states (DOS) function of the 
broadened Landau level n. 

The LL broadening results from impurity, phonon and 
electron-electron scatterings. In the experimental ultr- 
aclean GaAs-based 2D systems having mobility higher 
than lO'^m^/Vs the dominant elastic scatterings come 
mainly from residual impurities or defects in the back- 
ground rather than from remote donors and phonon 
and electron-electron scatterings are generally also not 
long-ranged because of the screening. On the other hand, 
since the microwave-induced magnetoresistance oscilla- 
tions occur at low temperatures and low magnetic fields 
in high electron density samples, the cyclotron radius of 
electrons involving in transport is generally much larger 
than the correlation length or the range of the dominant 
scattering potentials. In this case, the broadening of the 
LLs is expected to be a Gaussian form [e„ — {n + ^)uJc 
is the center of the nth Landau level], 

ImG„(£) = -(27r)5r-i exp[-2(e - Enf/T^] (9) 

with a B^/^-dependent half width expressed as 

r = (2c^,/7rT3)i/2, (10) 

in which t^, the single-particle lifetime or quantum scat- 
tering time in the zero magnetic field, depends on impu- 
rity, phonon and electron-electron scatterings. 

Expressions ^ and pH)) for the DOS of the nth level 
will be used in both the separated and overlapped LL 
regimes. The total DOS of a 2D system of unit area in 
the magnetic field, g{e) = — I'^^n(^)/''''^^B' '^^^ ^e 
written as the sum of the DOS without magnetic field, 
go, and a magnetic-field induced oscillatory part. When 
the LL widths 2r are much larger than the level spacing 
LUc (heavily overlapped LLs), the oscillatory part is much 
smaller than go and one can keep only its fundamental 
harmonic component that 

5(e)«go[l-25cos(27r£K)], (11) 

with S = exp(-7r2rV2w2) exp{-n /uJcTs)^ 

A reliable evaluation of single-particle lifetime Tj (or 
the LL width F), which must include effects of impu- 
rity, phonon and electron-electron scatterings, has not 
yet been available. We treat it as an empirical parame- 
ter, or, relate it to the zero-field linear mobility fiQ or the 
transport relaxation time th- using an empirical parame- 
ter a by 1/ts — Aa/rt^ — Aae/mfiQ. 

III. NONLINEAR RESISTIVITY AND 
MULTIPHOTON PROCESSES 

These formulations are convenient for current-driven 
nonlinear magnetotransport of any configuration. For an 



isotropic system where the frictional force F is in the 
opposite direction of the drift velocity v and the magni- 
tudes of both the frictional force and the energy-loss rate 
depend only on v = \v\, we can write F{v) — F{v)v/v 
and W{v) — W{v). In the Hall configuration with veloc- 
ity V in the x direction v =^ (u, 0, 0) or the current density 
Jx = J = Ngcv and Jy = 0, Eq. Q yields the transverse 
resistivity Ryx = B/Nge, and the longitudinal resistivity 
Rxx and differential resistivity r^x at given v as 

Rxx = -F{v)/{NU\), (12) 
Txx = ~{dF{v)/dv)/(N^,e^). (13) 

In the linear case, the resistivity Txx oscillations (peak- 
valley pairs) result from the resonant real photon-assisted 
electron transitions between different Landau levelsi^ 
We denote a real-photon assisted process in which an 
electron jumps across I LL spacings with the assistance 
(emission or absorption) of n photons as nLo:lu)c, or n:l. 
This process contributes, in the r^x^^/^c curve, a pair 
structure consisting of a maximum and a minimum on 
both sides of lo/lOc = l/n. 

Thus, the single-photon process 1:1, the two-photon 
process 2:2, the three-photon process 3:3,... all contribute 
to the maximum-minimum pair around uj/u>c = 1; the 
two-photon process 2:1, the four-photon process 4:2,... 
all contribute to the maximum-minimum pair around 
lo/ljJc = 1/2; the three-photon process 3:1, the six-photon 
process 6:2,... all contribute to the maximum- minimum 
pair around uj/oJc = 1/3; etc. The virtual photon pro- 
cess (zero-photon process) during which the electron may 
absorb (emit) one or several photons and then emits (ab- 
sorbs) them before finishing an intra-LL transition (as- 
suming separated LLs for convenience), which can be de- 
noted as 0:0, gives rise to an overall resistivity suppres- 
sion without structure.'^^ 

The effects of a finite current show up as the energy 
(frequency) shift lo^ = q\\ ■ v in the density-correlation 
function //2(Q|| , ^) in Eq- ©. This energy shift indicates 
that an electron, having momentum q|| and participat- 
ing in the system integrative drift motion of velocity 
v, carries an extra energy q\\ ■ v during its transition. 
In the case of low temperature (Tg much less than the 
Fermi energy ep) and large Landau- level filling factor 
{ly = Epji^c ^ 1), the major contributions to the sum- 
mation in Eq. ^ come from terms n ^ n' v, then the 
function C„.„' (x) has a sharp principal maximum near 
X = A:V. Therefore, as a function of the in-plane mo- 
mentum 1/2(911, ^) sharply peaks around g|| ~ 2Kp. 
As the composite effect of this energy shift, the electron 
possesses an extra energy 

ujj = 2Kfv ^ y/Sir/N.J/e (14) 

for its transition in addition to those provided by pho- 
tons. 

Therefore, the resonant condition applied in the linear 
case for electron scattering across / LL spacings with the 
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N=3.7x10'W^ =1200 cm^/Vs, a=4 




FIG. 1: (Color online) Magnetoresistivity Vxx vs ej — ujj/ujc 
at Lo/ujc = 1, 0.995, 0.99, 0.5, 0.49 and 0.475. 



assistance (emission or absorption) of n photons, nuj = 
lujc, should change to 

iUj ±nuj = ±luJc (n = 0, 1, 2, / = 0, 1, 2, ...) (15) 

when the system has an integrative motion of velocity 
V, i.e., a finite current density J — NsCv. We can still 
use symbol n:l to denote this process. We will also use 
symbol n- to represent all the 7i-photon participating pro- 
cesses (electron may jump across any LL spacings). 

Increase in the current density from zero reduces the 
rates of the scattering processes resonant in the linear 
case, such as those of zero-photon processes 0:0 every- 
where, of single-photon process 1:1 and two-photon pro- 
cess 2:2 at oj/uJc = 1, of two-photon processes 2:1 and 
four-photon processes 4:2 at uj/uJc = 1/2. But it may 
increase the rates of the scattering processes originally 
absent or very small in the linear case, such as those of 
single-photon process 1:1 and three-photon process 3:1 
at Lo/ojc — 1/2. 

At cyclotron resonance and its harmonics uj/uJc = 
1,2,4, the resonant condition ([15]) can be satisfied only 
at 

ej EE Ljj/tJc = 0, 1,2, (16) 

where the scattering rate peaks and the differential resis- 
tivity exhibits maximum. Therefore, at these magnetic 
fields the r^x oscillation with changing ej always has a 
periodicity Ae-,- « 1. 

The situation changes near subharmonic cyclotron res- 
onant fields. For instance, at oj/uJc = 1/2, the res- 
onant condition (|15p for zero-photon, two-photon, and 
four-photon processes (n = 0, 2 and 4) is satisfied at 
ej = 0,1,2,..., indicating these even-photon processes 
contribute resistivity maxima there. At ej = 0.5, 1.5, 
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FIG. 2: (Color online) Magnetoresistivity Txx vs ej at lj/uJc = 
1/3 and 0.5. Curves with 0-, 1-, 2-, 3- and m- are separated 
contributions from zero-photon, single-photon, two-photon, 
three-photon and multiphoton processes. 



on the other hand, the resonant condition is satis- 
fied for single-photon and three-photon processes (n — 1 
and 3), indicating these odd-photon processes contribute 
resistivity maxima. Because of this, the period of r^x 
oscillation with ej may shrink to Aej f» 0.5, or the oscil- 
lation frequency doubles, at lu/lOc = 1/2. 

At the third subharmonic field, lu/lOc — 1/3, in ad- 
dition to the appearance of resonant scattering of zero- 
photon and three-photon processes at ej = 0, 1 and 2, 
the resonance scattering can also appear at ej = 1/3, 4/3 
and 5/3 for two-photon processes and at ej = 2/3,4/3 
and 5/3 for single-photon processes, leading to possible 
shrinkage of the rxx oscillation period to Ae^ ~ 1/3, or 
the frequency tripling. 

These analyses are confirmed by the numerical results 
shown in Figs. 1 and 2, where we plot the differential re- 
sistivity rxx calculated from Eq. p3p as a function of the 
normalized current density ej — coj/ujc for a GaAs-based 
system of Ns = 3.7 x lO^^m'^ and = 1200mVVs 
at r = 1.5 K, irradiated by a 27 GHz microwave having 
incident electric field amplitude Eis = 3.8V/cm. The 
elastic scatterings are due to a mixture of short-range 
and background impurities, and the broadening param- 
eter a = 4, or the LL width 2r/uJc 0.72(w/wc)^/^, in 
the separated LL regime for uj/lUc < 1. The predicted 
rxx behaviors near u/lOc — 1 and lu/lOc — 0.5 are in 
excellent agreement with recent experimental findings.— 
The Txx oscillations at lu/lOc = 1/3 exhibit the observ- 
able frequency tripling under this strength of microwave 
irradiation. 
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IV. SUMMARY 

In summary, the electron in a dc-driven system having 
an integrative drift velocity v carries an extra average en- 
ergy ojj = 2Kpv in its transition between LLs, giving rise 
to a current-dependent resonance condition. The peaks 
of the differential resistivity r ,.,, arise from resonant mul- 
tiphoton assisted scatterings of electron jumping across I 
LL spacings. The current-induced consecutive emergence 
of different resonant multi-photon processes results in the 



frequency doubling and tripling of the resistance oscilla- 
tions near the 2nd and 3rd subharmonic fields. 
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